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Abstract 

We present a lattice calculation of the hadronic vacuum polarization and the lowest-order 
\ hadronic contribution (HLO) to the muon anomalous magnetic moment, = (g — 2)/2, using 



2 + 1 flavors of improved staggered fermions. A precise fit to the low-g 2 region of the vacuum 
polarization is necessary to accurately extract the muon g—2. To obtain this fit, we use staggered 
chiral perturbation theory, including a model to incorporate the vector particles as resonances, 
and compare these to polynomial fits to the lattice data. We discuss the fit results and associ- 



ated systematic uncertainties, paying particular attention to the relative contributions of the pions 
and vector mesons. Using a single lattice spacing ensemble generated by the MILC Collaboration 
(a = 0.086 fm), light quark masses as small as roughly one-tenth the strange quark mass, and 
volumes as large as (3.4 fm) 3 , we find a^ LO = (713 ± 15) x 10~ 10 and (748 ± 21) x 10" 10 where 



the error is statistical only and the two values correspond to linear and quadratic extrapolations in 
the light quark mass, respectively. Considering various systematic uncertainties not eliminated in 
this study (including a model of vector resonances used to fit the lattice data and the omission of 
disconnected quark contractions in the vector- vector correlation function), we view this as agree- 
ment with the current best calculations using the experimental cross section for e + e~ annihilation 
to hadrons, 692.4 (5.9) (2.4) xlO -10 , and including the experimental decay rate of the tau lepton 
to hadrons, 711.0 (5.0) (0.8)(2.8)xl0" 10 . We discuss several ways to improve the current lattice 
calculation. 

PACS numbers: 12.38.Gc, 12.39.Fe, 13.40.Em, 14.60.Ef, 14.65.Bt 
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I. INTRODUCTION 



The current Standard Model of particle physics has had tremendous success in describing 
fundamental properties of Nature. Currently, precise determinations of physical quantities 
need to be performed in order to find discrepancies, if any, between theoretical predictions 
and experimental measurements. Many, but not all, quantities agree quite well between 
theory and experiment. In those cases where there are discrepancies, there is a need to 
reduce the errors on theoretical and experimental determinations of said quantities to see if 
these discrepancies are significant. 

An excellent candidate for this improvement is the anomalous magnetic moment of the 
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muon, = (g—2)/2. The experimental world average is a° xp = 116 592 080(5.4) (3.3) x 10 

The incredible precision of this measurement arises from the recent experiments 
at Brookhaven National Laboratory This number deviates from current theoretical 
predictions by 0.7 to 2.7 standard deviations, depending on the method of obtaining the 
Standard Model number (see the recent reviews 3|, |4j). 

The largest contribution to the error of the theoretical calculation lies in terms arising 
from QCD. The electroweak sector can be reliably calculated using perturbation theory but 
this is not the case for QCD. The first QCD contribution comes in at order 0(a 2 ) (a = e 2 /47r 
is the fine structure constant) shown in Figure [1] which is the hadronic contribution to the 
photon vacuum polarization. The next QCD contributions come in at C(a 3 ), the light- 



by-light contribution, which is currently being studied using 



calcu 
Ref. 



ations (e 



■9-i R Q, [§], Isj]), and chiral perturbation theory h| (for a recent review, see 



attice techniques {5], model 



111]). We focus here on the larger of these two contributions, that coming from the 
QCD corrections to the photon vacuum polarization. 

Since perturbation theory is unreliable for such hadronic contributions, for a long time 
one has relied on a dispersion relation to relate the imaginary part of the photon vacuum 
polarization to the real part [k]] which is needed in the diagram shown in Figure 
using the experimental e + e~ —* hadrons total cross-section to obtain the imaginary part. 
Additionally, one can relate the decay r — > hadrons to the e + e~ cross section in the isospin 
limit 3|, |4| for a part of the energy range. Whether or not one includes the r data gives the 
range of deviations discussed above between theory and experiment. 

A purely theoretical calculation of a^ LO is desirable for at least three reasons: first, it 
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is very important to corroborate the dispersive/experimental analysis with an alternative 
calculation in order to have confidence in the Standard Model comparison with experiment, 
and second, these calculations may eventually compete with or surpass the precision and 
accuracy of the dispersion relation/experiment combination. Finally, from a theoretical 
perspective it is quite interesting. Such theoretical calculations may in the future provide 
even more stringent tests between the Standard Model and experiment, a desirable outcome 
in the LHC era. 

At present the only way to compute aj^ LO theoretically in a model-independent way is 
to apply lattice gauge theory techniques to calculate the hadronic vacuum polarization. 
This has already been done in the quenched approximation with domain wall fermions [jjj] 



and improved Wilson 
staggered fermions 



16. 



ermions 



151 ]. as well as with dynamical (2 + 1 light flavors) "asqtad" 



171 ]. The main difficulty is that the largest uncertainties occur in 



the low-momentum region which gives the dominant contribution to a^ L °. The goal of the 



current work is to incorporate staggered chiral perturbation theory (SXPT) 18[ coupled 
to photons to derive a physically based fitting function to accurately describe this low- 
momentum region. Initial results from this fitting technique were presented in 19] , and to 
our knowledge are the first chiral perturbation theory results for the vacuum polarization 
presented in the literature. 

At this stage we must point out that an expansive and lively debate has recently taken 
place in the literature, and the lattice QCD community in general, concerning the theoretical 
validity of the staggered fermion formulation when applied to numbers of quarks different 
from a multiple of four, as in the present case of 2+1 flavors. We refer the interested reader 



to the most recent papers 



Lattice 2006 



21 



22 



23 



24 



251 ] . and also to the review by Sharpe at 



261 ] . The important point stressed in the review is that staggered fermions 



are non-local for non-zero lattice spacing, but that the non-locality is a lattice artifact that 
vanishes in the continuum limit, provided the quark mass is not taken to zero first. Since 
we work at fixed lattice spacing, this may be cause for concern. As we shall see, the results 
presented below do not indicate signs of non-locality, and our future plans include taking 
the continuum limit when ensembles with smaller lattice spacing become available 1 . We also 
plan to repeat our calculation with 2+1 flavors of domain wall fermions as a final check. 



1 The MILC collaboration is currently generating such an ensemble. 
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This paper is outlined as follows. We discuss in Sec. [TT] the basic details of the lattice 
calculation of the vacuum polarization and how to use this to extract the lowest order 
hadronic contribution to the muon g—2. In Sec. IHII we discuss staggered chiral perturbation 
theory coupled to photons as external fields, and compute the one-loop pion and kaon 
contribution to the photon vacuum polarization. We find here that this is not an adequate 
description of the lattice data, and in Sec. IIVI we incorporate into the chiral perturbation 
theory the lightest vector fields, whose contribution to the vacuum polarization we find 
dominates over the pions. Sections [V] and [VI] include our discussion of the fits, final results 
for the muon g — 2 for this lattice spacing, and systematic uncertainties in the calculation. 
We conclude our study in Sec. IVIH 



II. LATTICE CALCULATION OF THE MUON g-2 

In the continuum, Euclidean and gauge invariance along with the Ward-Takahashi iden- 
tity allow us to write the vacuum polarization tensor in momentum space as 

IV(g) = (q\u - q,q,mq 2 ) • (1) 
We can use II(g 2 ) as defined above in the space-like region to calculate the lowest order 



contribution to the muon g — 2 from the expressions 



a 
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oo 



C° = F dK 2 f(K 2 )U(K 2 ) , (2) 



o 



m 2 K 2 Z 3 (1 - K 2 Z) 
f( R ) = -i I L 272 , (3) 

K 2 - JK 4 + AmiK 2 

2m 2 K 2 ' K ' 

fi(K 2 ) = ^Y^n^-n^o)] . (5) 

i 

Here Qi is the electric charge in units of e, and IT^i^ 2 ) the vacuum polarization for the 
i th quark flavor with space-like momentum-squared K 2 . U(K 2 ) is logarithmically divergent 
and is renormalized by subtracting its value at K 2 = 0. Since f{K 2 ) diverges as K 2 — > 0, 
we see immediately that the largest contribution to the muon g — 2 will come from the 
low momentum region; in particular the dominant contribution to the integral comes from 
momenta near the muon mass, K 2 w m 2 . 
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The vacuum polarization for the i th quark flavor in the continuum is given by (no sum 
on i) 

K(D = J d'xe^-y^J^JKy)) , (6) 

where J % A%) = ipi(x) , ~f^ipi(x) is the vector current. 

For staggered quarks, the conserved current has the following form, in terms of the one- 
component staggered fields x, 

= 7^0) [Xii x + P)Ul(x)x(x) + Xi{x)U„{x)x{x + p,)] , (7) 

where t]^{x) = (— \y-"^<y. x " is the staggered fermion phase, and satisfies a discrete conserva- 
tion law 

= (8) 

(up to possible contact terms depending on the correlation function) where A M is a backward 
difference operator. This implies an exact lattice Ward-Takahashi identity for the two-point 
function, given by Eq. ([I]) with — > q^ = 2sin(ag At /2)/a. It is important to note that 
this Ward-Takahashi identity must be satisfied (in a numerical simulation) on each gauge 
configuration separately (to numerical precision). This is true for all of the calculations 
presented here and serves as a strong check of the calculation. 

In evaluating the quark propagators using the current above, we do not include the Naik 



2a 



29 



30| used by the MILC 



term [27j appearing in the improved staggered (asqtad) action 
collaboration to generate the gauge field configurations used in this study. This is because 
the form of the Ward-Takahashi identity in the presence of the Naik term is not as simple 
as it is for the nearest neighbor point-split current: an extra term appears. The effects of 



not including the Naik term were discussed in [la]. Recall that the Naik term is meant 
to improve the dispersion relation of the naive staggered fermion which breaks down for 
large values of the momentum. Thus this omission should not significantl y e ffect the long- 



distance, low momentum physics with which we are mainly concerned. In [3l| it was shown 
that the effect on hadron masses, and the vector and pion meson masses in particular, from 
the Naik term was negligible compared to the statistical error for the types of parameter 
values used in this study, i.e., small quark mass and lattice spacing. While we know of no 
significant effect due to this omission, it is possible that an unknown error is introduced by 



the "mixed action" we are using. Since 



161 ] showed little difference (within statistical errors) 



whether or not the Naik term is present, we assume for now that any error introduced by 
this mixed action is negligible. 

Further details of the lattice calculation can be found in , Q, Q, although we 

make one last comment here as it pertains to our discussion later. There are two types of 
Wick contractions (shown in Figure [2]) that need to be evaluated for a complete calculation; 
however the disconnected diagrams shown in Figure [2] have not been computed due to 
the well-known practical difficulty in calculating them. The disconnected diagrams, when 
summed over u, d, and s quarks, cancel in the SU(3) flavor limit, and they are suppressed 
according to Zweig's rule. As such we do not expect them to make a significant contribution 
to the final result, so for now we leave the calculation of the disconnected diagrams for a 
future work. 

Table [J shows the lattice parameters used in this calculation. We employ the "fine" 
MILC configurations (a = 0.086 fm) generated with the asqtad fermion action with the 
light valence quark mass set equal to the degenerate light sea quark mass and fixed value 
for the strange valence quark mass equal to the strange sea quark mass 2 . We will discuss 
the lattice data and the fits in Sec. [V] 

III. STAGGERED CHIRAL PERTURBATION THEORY WITH PHOTONS 

Since the low-energy, small quark mass regime is the most important part of the lattice 
calculation (because of the kernel in Eq. ([3])), we will need a precise, theoretically sound, 
fitting function to model the lattice data. In principle, if the lattice volume were large enough 
to have arbitrarily small momenta, and the quark masses took on their physical values, then 
any simple, smooth parameterization of U(K 2 ) would do, provided the statistical errors 
on the lattice data were small enough. Even on the state-of-the-art lattices used in this 
study, this is not the case, so the lattice calculation must be leveraged as much as possible 
with a good fit function. Such a function for the vacuum polarization can be determined 
in chiral perturbation theory (XPT) coupled to photons treated as external fields, as in 



34j . for example. Coupling XPT to photons is simple: to treat the photons as external 



2 This is nominally the "physical" value of the strange quark mass, although it was noted in [32 ] that this 
value is slightly too high on these lattices, and intstead am s sa 0.027 corresponds to the physical strange 



quark mass. This difference should, as in Refs. 32j, [33j, have a negligible effect 
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fields we merely replace the ordinary derivative with a covariant one. The generalization to 
staggered XPT (SXPT) is also relatively straightforward, but there may be additional lattice 
artifacts ("taste-symmetry 3 violations") that arise from the interactions of the pions with 
the photons. As we shall see, these terms do not arise at the order to which we work. 

First let us write down the staggered chiral Lagrangian in the presence of photons. The 
chiral symmetry for three flavors of staggered quarks is an 577(12) £ x SU(12)r symmetry. 
We incorporate the pseudo-Goldstone bosons in the field E, where E — > LHW under a chiral 
transformation, with L e <S77(12)£ and R G SU (12) r. 

In Minkowski space, we have the Lagrangian 

f 2 u f 2 

C 2 = ^-Tr^E^E)*] + ^-Tr[.ME + T) M] - a 2 V , (9) 
8 4 

where V is the taste-symmetry breaking potential for multiple flavors of staggered quarks 
jisj ]. M. = diag(m u l4x4, m^l^, m s i4x4) is the light quark mass matrix, / is the tree-level pion 
decay constant (normalized here so that the physical f n ~ 131 MeV), and /i is a constant 
with dimensions of mass. The subscript on £ 2 signifies the 0(p 2 , m 2 , a 2 ) chiral Lagrangian 4 . 



The explicit form of V is given in [18[, and for our purposes gives rise to the splitting of 
the various taste mesons. It additionally gives rise to two-point interaction terms which mix 
the flavor-neutral vector- and axial-taste mesons, but these will not have any effect on the 
present calculation. The covariant derivative is defined as 

D tt E = d tt E + ieA tt [Q 1 H\, (10) 

where Q is the light quark charge matrix in units of e, and is the photon field. 

Note that we have dropped the term in the Lagrangian which is proportional to m 2 ,. 
This term, due to the anomaly, gives the taste- and flavor-singlet meson (rjj) a large mass 
which we send to infinity, as generally the effects of the rl l are decoupled from the low energy 
regime. In addition, the r)j, along with all the other flavor-neutral pions, does not contribute 
to the photon vacuum polarization. 

The question now is whether or not there will be additional terms which can violate the 
taste symmetry at the order to which we work, but which include photon interactions. Due 



3 Each staggered fermion corresponds to four degenerate Dirac fermions in the continuum limit. These are 
called tastes, in analogy to flavors. 

4 Generally, for staggered fermions the leading discretization errors begin at 0(a 2 ). 
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to the gauge structure of the theory, the only place where the gauge fields can enter is either 
through a derivative term, or by inserting terms that contain the field strength tensor F^ v 
in a Lorentz- and gauge- invariant way. The first thing to note is that V does not include 
derivative interactions, so there are no photon terms arising from covariant derivatives in V. 
Terms with only a single factor of F^ u are not allowed due to Lorentz invariance. Thus, in 
the 0(a 2 ) La grangian, there are no taste-violating terms that contain photons. 

As usual in chiral perturbation theory, working at one-loop order in £2 requires us to 
include tree-level terms from £4, as these two sets of terms are of the same order in power 
counting. Additionally, these tree-level terms absorb the infinities (and scale dependence) 
that arise from the one-loop integrals. The only terms from £ 4 that violate the taste sym- 
metry are 0{a 2 m 2 K} _^'p 2 , a 4 ) since we must include factors of a 2 . Terms in the continuum 



0(p 4 ) Lagrangian 34| also contribute at this order. 



In 35| all taste-symmetry breaking operators at 0(a 2 m 2 , a 2 p 2 , a 4 ) which include external 
fields along with the meson field have been enumerated, and thus should include all possible 
terms that can contribute in the present case. It is easy to see that none of the terms listed 



in 35[ contributes here, just as many of the 0{p 4 ) terms do not. The reason is that any 
terms which include a covariant derivative on E vanish when E = 1 in the 0(a 2 m 2 , a 2 p 2 , a 4 ) 
terms 5 (there can not be any factors of momenta if each term including a field strength 



tensor is to have the correct dimensions). Thus, all operators listed in 35j vanish, and from 



the 0(p 4 ) terms 34j, we only have the following terms which will contribute 

£4 = LioTrp^EFH + 2# 1 Tr[F Mi ,FH . (11) 

We can easily see that both of these terms will contribute in the same fashion when we set 
E = 1. 



Terms with only external fields (contact terms) were not considered in 35], although in 
the current calculation they could play a role, just as the Hi term in Eq. (flT|) does. We 
can see that explicit taste-symmetry breaking terms of this form cannot arise. Any terms 
included here have to be 0(a 4 ,a 2 p 2 ), but two factors of F^ u mean any such contact term 



5 Recall that we are calculating the photon two-point function, so to obtain a tree-level term which will 
contribute, then we can only have factors of the photon field. Thus, we will ultimately set £ = 1 in all 
£4 terms. 
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begins at 0(p 4 ) 6 . Thus, at most one power of can appear in each term. Lorentz-invariant 
terms with one factor of are not allowed, as they require more derivatives. However, it 
could be possible to break the Lorentz symmetry by coupling the Lorentz indices on the field 



35 



Following 



strength tensor to the vector indices on the taste matrices, as in Refs. 
a similar spurion analysis, we find that the only terms that would arise are those which 
involve the diagonal elements of the field strength tensor, which all vanish. Thus there are 
no 0(a 4 , a 2 p 2 ) taste violating terms. We could get 0(a 2 ) corrections to various parameters 
(such as L 10 or Hi), but these would not arise from violations of the taste symmetry, but 
instead from other lattice discretization errors. 

While the above argument is purely technical in nature (enumerating all terms that are 
allowed to appear in the theory), we can see from a more physically intuitive picture that 
there are no taste violations (at least to this order) from photon interactions. Taste violations 
generally arise from gluonic exchange with momentum near the cutoff. These could also arise 
from photon exchange with momentum near the cutoff, but would be suppressed compared 
to the gluon exchange by a factor a/a s . More importantly, however, in this lattice calculation 
these interactions do not occur because we do not treat the photons as dynamical fields. 
In other words, our simulation is a quenched QED simulation, and as such we cannot have 
such taste violations arise. This does not mean there are not any taste violations here; they 
can, and do, arise in the meson mass splittings which we already know come only from high 
momentum gluon exchange. 



A. One-loop calculation with pions 



Using Eq. (Q, it is straightforward to calculate the one- loop contribution to the photon 
vacuum polarization. As usual we write £ = exp [i <£//], where $ is the 12 x 12 matrix of 
pseudo-Goldstone bosons 



u 


7T + 


K+ 




D 


K° 


K~ 


K° 


S 



(12) 



Note the external photon field is 0(p) in the chiral power counting scheme 34 ] . so ~ d^Ay — d v A 
is 0{p 2 ). 
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where $ is written in the flavor basis so U is a uu meson, D is a dd meson, and S is an 
ss meson. Also, each element shown in Eq. f|T2|) is a 4x4 matrix, written as, for example, 
U = Ea U aT a , where 

T a G {{ 5 , z£ m5 , z^, f M , 14x4} (13) 

is a Hermitian taste matrix (a = 1 . . . 16), 1^4 is the 4x4 identity, the £^ are Hermitian 
Dirac matrices which satisfy {C^,Cv} = 25 M1/ , ^ v is defined as (similarly with £ m5 ), and 

< ^ in Eq. ffl~3"l) . Note that although we are in Minkowski spacetime, we use Euclidean 
taste matrices as in 36^ . 

Expanding E to quadratic order in $ gives two types of couplings to the photon, as is the 
usual case for scalar QCD. Both vertices have two pions, and one vertex has a single photon 
while the other has two photons. From these vertices, we see that there are two possible 
diagrams, shown in Figure [31 where the mesons in the loops are only the charged pion or 
kaon, of any of the 16 tastes. Note that although the second diagram is a constant in q 2 , 
where q is the external momentum, it is required for gauge invariance of the total amplitude 
at this order. 

The calculation for these two diagrams is rather straightforward, and we get the following 
result, after continuing the external momentum to Euclidean space, q 2 — > — q 2 , 

,2 \ o 1 / m 2 



M,t 

where 



(14) 



^-^Kttt^t)-! ■ (15) 

M e {ir + ,K + }, and the sum over t is over the 16 tastes of mesons. A is the chiral cutoff 
scale. The subscript 4 refers to the fact that we have not yet taken into account the reduction 
from 4 — >• 1 tastes per flavor which must be done to arrive at the physical number of mesons. 
The mass for a meson M t consisting of two quarks a and b and with taste t in terms of the 
Lagrangian parameters is 

m 2 Mt = fi(m a + m b ) + a 2 A t (16) 
where At is a combination of various parameters from V, and can be found in 18[. The 



measured values on the lattice can be found in 321 ] . 

The above expression is the lowest order result for the vacuum polarization in chiral 
perturbation theory. In particular, to this order, it is a valid description of QCD and needs 
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no alterations whether the disconnected diagrams are included in the lattice calculation or 
not. Imagine a lattice calculation in the SU(3) flavor limit; the disconnected diagrams in 
the electromagnetic current correlation function cancel exactly. However, in Eq. Othe only 
change is to set = tuk- In other words, effects from these diagrams must vanish like 
(m£ — m s ) 2 and hence are higher order. 

So far, Eq. ({141) includes the effects of each flavor of quark having four tastes. In the 
simulations, the fourth root of each quark determinant is (formally) taken to remove the 

nn 

extra unphysical tastes, so we need to include these effects by hand in SXPT [18|, 137J . To do 



this, we use standard quark flow techniques (for example, [18] and references therein), with 
minor adaptations due to the presence of the external photon field. 

As mentioned before, the interaction terms involve two pion legs and either one or two 
photon legs. The quark-level diagrams for these two vertices are shown in Figures H] and O 
Combining these vertices in such a way to correspond to the diagrams in Figure [3j we see 
that both diagrams include two internal quark loops. Thus, every term in our result must be 
multiplied by a factor of (1/4) 2 = 1/16 to account for the transition from 4 — >• 1 tastes per 
flavor. 7 Incorporating the necessary factors of 1/4 in our expression above, and correcting 
for the exclusion of the disconnected diagrams, we have 



n l -loo P(g , 




G( 4^\_8 + l ln ^4 t 



q 2 J 9 3 \ A 2 



A . (17) 



The constant A includes the analytic terms from the 0(p 4 ) Lagrangian (shown in Eq. (Ilip ) 
as well as terms that can arise from non-taste- violating 0{a 2 p 2 ) operators mentioned above. 
Since the current work is only done at a single lattice spacing, we do not worry about the 
explicit a 2 dependence in A, although calculations at multiple lattice spacings would require 
writing 

A = Aq + Axo 2 , (18) 



7 Note that for the vacuum polarization computed in the chiral theory all quarks involved are "sea" quarks 
since we compute the photon two-point function, even though in the lattice calculation there is a distinction 
between sea and valence quarks. Thus, although one might assume that some quark lines are treated as 
valence quarks and do not get a factor of 1/4, this is in fact not the case. This is also the reason why, 
although we exclude disconnected diagrams in the lattice calculation, wc include what look like, at the 
quark-flow level, disconnected diagrams. Both of these statements can be made concrete by treating the 
XPT calculation in a partially quenched setting, and one finds the same results presented here. 
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with Aq oc Lio + 2Hi, and A\ some coefficient parametrizing the discretization errors. Note 
that this additional discretization error does not affect the muon g—2, as it will be subtracted 
off when renormalizing the final result. The only relevant discretization errors come from 
the different meson masses in Eq. (1171) . 

We would like to examine Eq. (TTTf) in both the zero momentum and the chiral limits. 
These correspond respectively to x — ► oo and x — > in the function G(x). The first limit is 
well defined, and we have that 

lim G(x) = - . (19) 

x— >oo 9 

However, in the chiral limit, when x — ► 0, the function G(x) diverges like 

\imG(x) = -In ( - J +0(x) . (20) 
x^o 3 \x J 

This divergence is exactly canceled in the full vacuum polarization by the term | ln(m|^ t /A 2 ), 
so n}~ loop (g 2 ) is finite in the chiral limit, as we expect. 

We find immediately that this expression does not describe the data at all. We can fix 
the overall constant A by matching the data at the lowest value of the lattice momentum, 
q 2 = {2n/aT) 2 rs 0.0223 GeV 2 , for the largest lattice used in our study and use the values 



of the different meson masses from 



321 ] . After fixing the overall constant there are no 



free parameters, and the resulting fit function and the lattice data for light quark mass 
ami = 0.0031 are shown in Figure [6j The slope for low q 2 is smaller by roughly two orders 
of magnitude in the fit function compared to the data. It is obvious that this cannot be a 
complete description of the low-energy physics for the photon vacuum polarization. 

We can use our knowledge of the quenched simulations Q, [jjj] to understand this. In 
the quenched case, the entire result for the vacuum polarization is described by a form 
dictated by vector particle intermediate states since they are absolutely stable in this case. 
We imagine that the vectors can also make a significant contribution to the photon vacuum 
polarization in the dynamical case, so we should include them in our calculation. Thus, the 
next section discusses how we incorporate vectors in SXPT. 



IV. INCLUDING VECTORS IN SXPT AS RESONANCES 



Including the lightest vector mesons such as the p or the u is more tricky than includ- 
ing the light pseudoscalars. Unlike the pions, the light vectors have no special role when 
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discussing chiral symmetry. This is not an inherent problem, although there could be some 
model dependence which arises depending on how one includes the effects of the vectors. 
Here we use the formalism of Ecker, et al. in 38j, treating the vectors as resonances in the 



effective field theory. This describes the vectors with antisymmetric Lorentz tensor fields, 
and not with Lorentz vector fields. With the correct free-field Lagrangian, this will enforce 
the correct number of degrees of freedom, and is useful when coupling the vectors to the 



photon 



The particular form of the vector field should not alter the physics, as was 



discussed in Ref. 



The octet of light vector mesons is an antisymmetric spin-1 field, V^,, which is an SU(12) 
matrix. The continuum free-field Lagrangian in Minkowski space is 



C v = — Tr 
2 



d^ v d p Vf> v - \m 2 v V^V^ 



(21) 



where my is the vector meson mass. 

From the point of view of quark content, the vector fields have the same status as the 
light pseudoscalar fields. For example, a p + has the same light quark content as a 7r + , 
although the relative spins of the constituent quarks are different. There is a fundamental 
difference from the point of view of chiral symmetry, in that the vectors do not have any 
special properties when the light quark masses vanish. What is important, however, is that 
the vectors and the pions have the same transformation property under the unbroken vector 
SU (12) subgroup, as this is governed solely by the quark content. The primary difference 
is that we cannot determine the leading dependence of the vector mass my on the quark 
masses from chiral symmetry alone, as we can for the pions. 

As is standard in XPT, we expect the SU(12)l x SU(12)r chiral symmetry to sponta- 
neously break down to the vector subgroup SU(Yl)y. It is under this subgroup that the 
vector and pseudoscalar fields will have the same transformation property 

s -> u v m v , 

V, u - UyV^Uy , (22) 



By including the vector resonances, we emphasize that we are now resorting to a model, so in principle 
there can be model dependence in the results. However, at present this model is required to adequately 
fit the low momentum region of the lattice data and as such need only furnish a smooth function that 
accurately describes that data. When larger volumes are used, one might dispense with such fits as the 
results could be reliably extracted from the actual data if enough low-momentum points are reached. 
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where Uy G SU(12)y is a constant vector transformation. With regard to the underlying 
chiral symmetry, any transformation which reduces to Eq. ff22l) in the appropriate limit will 
lead to the same physics. A simple choice for the underlying transformations is made by 
defining a 2 = £, where under SU(12)l x SU(12)r we have 



a -> LaU(x) ] = U(x)aR ] 



(23) 



U(x) is an element of SU(12), and since it depends on the pion fields, we see that it is 
position-dependent (the explicit form of U(x) will not be needed). We can now define a 
transformation rule for the vector octet as 



(24) 



When SU (12) L xSU (12) R -> SU(\2) V , we set U(x) = L = R = U Vl and all the fields have 
the correct transformations under the unbroken vector subgroup. 
The interaction terms to lowest order are given by 



f\ 



iG 



v 



(25) 



2V2 " " UJ+ ' 2V2 
where fy and Gy are real couplings (fy is the tree-level vector decay constant). We use the 
definitions 



a' 



= eQ(d^A u - d v A») . 



(26) 
(27) 
(28) 



The second term in Eq. ([25]) will not affect the present calculation. An explicit form for the 
vector field is given by 

p+ K* + N 

P~ + K * 



( i n o 1 1 , , 



V, 



[US 



'*() 



V 



K* 



K 



*0 



(29) 



where in this expression, each entry is a 4x4 matrix, corresponding to the 16 tastes of each 
vector meson. 

We pause our discussion for a comment about the different tastes of a staggered vector 
meson. There are sixteen tastes, falling into eleven multiplets which are irreducible repre- 



sentations of the lattice symmetry group 



401 ] . Each of these irreps, in principle, will have 
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different masses, but empirically the mass differences of the tastes are negligible. We use 
this fact to ignore these mass differences, to a first approximation, in the current analysis. 
Of course this does not imply that there can be no significant effects from taste-changing 
interactions among the vector particles in general, although we assume they give rise to 
errors which are smaller than the statistical and systematic errors already present in our 
calculation. 

Finally, we should include terms in the vector part of the Lagrangian that could violate 
the taste symmetry, but we choose to neglect them for the following reasons. On one hand, 
taste violations arise through additional terms in the masses of the vectors (as happens with 
the pions), and these we neglect because empirically they are negligible. The other source of 
taste violations is additional terms that do not contribute to the tree-level masses but arise 
in one-loopcalculations just as the flavor-neutral hairpins arise in pion masses and decay 



constants 



41] . We expect these contributions to be similar in size to the mass splittings 



(as in the pion sector), and thus negligible. 
A. Tree-level vector contribution 

The tree-level calculation of the vacuum polarization involving the vectors is straight- 
forward. We set a = 1 in the above expressions and find the vertices which couple the 
vectors to the photon. There are two two-point vertices: one where the photon couples to 
the p° and another where the photon couples to the u s . The relevant diagram is shown in 
Figured and from it, we get the following contribution to the photon vacuum polarization, 
after continuing the external momentum to Euclidean space, 



a (4vr/ v 



47T 



+ 



p 2 + my p 2 + rriy 



A7r3p 2 + m 2 v v ; 



In this expression, we have not included an explicit taste index on the vector mesons since we 
are neglecting the splittings between the various tastes. The full analysis requires making 
the transition from four to one tastes per flavor. By examining the Lagrangian above a 
similar analysis to the pion case leads to an overall factor of 1/16 times a sum over the 16 
tastes of vector mesons. Neglecting taste violations in the vector masses, the 16 degenerate 
contributions cancel the factor of 1/16. 

If we were to measure the decay constant as well as the vector masses on the lattice, then 
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this expression has no free parameters. Notice that at this (lowest) order, the masses of both 
the p° and u)% are the same, my. We can set this equal to the p mass, and any discrepancies 
enter at 0{jP 
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A quick analysis shows that for small momentum, the tree-level p contribution dominates 
the one-loop pion contribution above. Expanding Eq. (1301) around p 2 = 0, we get 

nr (p 2 ) = nr (o) + ^^IyLiL + 0(p «) , (si) 

47r 3 mr p0 
while for the one-loop pion contribution we get 

n£"V) - n-'-(o) + £± Y. 30&7 + °^ ■ < 32 > 

M,t M * 

To more easily compare these two expressions, let us first subtract off the zero momentum 
contribution, because it is the slope of the low-momentum region that we are interested in. 
Taking the ratio of these gives 

Iltwe(p2) _ TI£ CC (0) 640(47r/V) 2 



rW) - 11,(0) r 



(33) 



Using a rough estimate for fy m 200MeV 42j], we see that only for very light pion masses will 
the pion contribution dominate. For the ensemble with the light quark mass ami = 0.0062, 
we see that 



200GeV" i , (34) 



and so for a vector mass of ~ 800MeV, Eq. (1331) gives about 50. In the practical simulations 
we expect the vector contribution to dominate. Even at the physical masses in the continuum 
limit, this ratio is still roughly 10, so although the pion contribution becomes more important 
for lighter quark masses, the majority of the slope of the vacuum polarization comes from 
the vectors. 

Although the contribution of Eq. (1171) is much smaller that of the p, we will include it 
to determine the effects it has on the final fits. As such, we will have to include also the 
one-loop diagrams which involve the p for consistency. This is done in the next section. 



B. One-loop terms with the p 



The only diagram that will contribute here is shown in Figure [8] (there is a second, with 
the pion tadpole attached to the other vertex). There is also the possibility of a term 

16 



arising as a self energy correction to the vector (so the vector line would have a pion loop 
inside). This has been evaluated (Eq. (28) of 43j), but it is 0(p 4 ). Although this term is 
proportional to an unknown low energy constant that in principle could be rather large, to 
include it consistently would require us to include higher order terms from the pure pion 
sector. As such, it is consistent to neglect this piece. 

A straightforward calculation gives the one-loop terms. Putting this together with the 
tree-level vector contribution we get the complete one-loop result for the vectors, 



n 



1— loop 



(p 2 



a (47r/ 



Vj 



4vr 3 
1 



+ 



p 2 + m v 
12 



4 



(4tt/) 2 16 



p 2 + mi 



where the usual chiral logarithm is 



(4vr/) 2 16 



2 N 



+ 



p 2 + m v 



■(3me + m s 



(35) 



£(m 2 



(777, 
~j\2 



(36) 



with A the chiral scale, as before. The factors of 1/16 arise in the same way as in the purely 
pion quark-flow analysis. As we are neglecting taste violations from the p, we can treat it 
as a taste singlet, and then the p-27r-photon vertex is identical in quark-level terms to the 
2 photon-2 pion vertex. Thus, the quark-level diagrams are identical to those in Figure 



with one of the photon lines taken to be a p line. We again keep my 



m n o 



Tfi u] , as 



. We have 



the discrepancies are 0(p ) and are still higher order than this expression 
included an unknown coefficient, Cy, which would arise from analytic terms from the higher- 
order Lagrangian. This is necessary to cancel the scale dependence coming from the chiral 
logarithms. 

We can see that this is roughly the same order of magnitude as the one-loop pion contri- 
bution, so if we were to include the one- loop pion expression in our final result, we would 
in principle need to include this. We do not explicitly include a constant term in Eq. (1331) . 
even though to this order it would be consistent to include such a term. Following the 
reasoning as in the pion sector, such a term would arise in the same manner. Since the 
full one-loop result would include Eqs. (1T7|) and (1331) . this would amount to a change in the 
definition of the constant A in terms of the underlying chiral parameters. As we are not 
explicitly writing down the higher-order terms in the Lagrangian in the p sector, nor are 
we trying to extract the values of specific low-energy constants, we only need the mass and 
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momentum dependence. Additionally, in the end this does not matter since to determine 
the muon g — 2, we subtract off the value of the vacuum polarization at zero momentum, 
so any overall constant does not affect the physics. Note also, though, that the one-loop 
corrections in Eq. (1551) could be absorbed into the definition of the parameter fy, as there 
is no momentum dependence in these terms. We will discuss this further in the next section 
when we discuss the fits to the lattice data. 



V. NUMERICAL RESULTS 

We plot the results for the 2+1 flavor vacuum polarization as a function of q 2 in Figure [H 
including a close-up of the low-g 2 region. As before HQ, I Il7j ]. for larger q 2 , it is independent 
of quark mass, but as q 2 — > 0, there is a significant difference in the three plots. IT(g 2 ) was 
computed using two quark propagator source times, and 48, and results on the same lattice 
were averaged. We did not compute with both sources on every lattice, however, and strange 
and light valence quark propagators were calculated on overlapping but different subsets of 
each 2+1 flavor lattice ensemble (see Table [I]). The statistical errors shown in Figure [H 
and quoted throughout, come from a single elimination jackknife procedure. Because the 
numbers of strange and light valence quark propagators differ within the same ami = 0.0031 
or 0.0062 ensemble, we adopted a jackknife procedure where each light or strange quark 
propagator calculation of IT(g 2 ) is treated as a separate measurement, instead of the 2+1 
flavor value of II(g 2 ). This is not expected to cause difficulty since the much heavier strange 
quark means the light and strange quark propagators computed on each lattice are roughly 
uncorrelated. Also, because of the electric charges, the light quark contribution is explicitly 
weighted five times more than the strange quark one. For the ami = 0.0124 ensemble, the 
jackknife error estimate for H(q 2 ) was calculated in the usual way. In all cases the errors 
were not significantly altered by increasing the jackknife block size to five, or by computing 
them using a simple binning procedure with bin size of ten configurations. 



In Figures [10] and [TT] we compare continuum three-loop perturbation theory 44| with the 
lattice calculation of — II(g 2 ) for < q 2 < 8 GeV 2 for mi = 0.0031 and 0.0124, respectively. 



The perturbative result is given in the MS scheme, and the bare quark mass has been 



matched using the renormalization factor given in 



451 ] . The results were forced to agree at 



/i = 2 GeV by a imposing a simple additive shift to the perturbation theory curve. The 
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lattice results for 3 < q 2 < 8 GeV 2 agree impressively with perturbation theory. For lower 
values of q 2 and mi = 0.0031, the lattice value increases faster until about 0.5 GeV 2 when 
the diverging perturbation theory result overtakes it again. For the heavier mass, the two 
results coincide until about 1 GeV 2 where the perturbative one again overtakes the non- 
perturbative one. This behavior with quark mass is mostly indicative of the variation with 
quark mass in the lattice results since the perturbation theory curve is insensitive to the light 
quark mass. This can also be seen from Figure [9j the two lattice results are clearly separated 
at 1 GeV 2 . The above comparison is, of course, not unique. We could have matched the 
results at a value other than 2 GeV. Still, it is clear from Figures [10] and dH that a large 
range of momentum values exist where perturbation theory is expected to be reliable, and 
the precise matching point will have little effect on the ultimate value of a M . The matching 
point of 2 GeV is a natural choice, given the lattice scale in our study. This comparison 
is a non-trivial check of the lattice calculation of the vacuum polarization; in particular it 
indicates the correct number of physical quark degrees of freedom are accurately accounted 
for when using staggered fermions. 

For the low energy regime, we have tried several different fits. In every case the fits are un- 
corrected but performed under the single elimination jackknife procedure described above. 
Following previous calculations of the muon g— 2 Q, H], we try various polynomial fits, 
up to quartic order in q 2 . The general form for these fits is 



with x = q 2 , and n ranging from 1 to 4. We will refer to these fits individually as "Fit 
n." For larger n these fits have adequate values of x 2 ) but tend to undershoot the lowest 
momentum points except for n = 4. For n < 3, the fits have very poor x 2 and do not fit the 
low-g 2 reqion well. In Figure [121 we display the cubic and quartic fits for each light quark 
mass. In Table Ull we show the results for the fitting parameters for cubic and quartic fits at 
each of the three different quark masses. Additionally, we use these results to determine the 
q 2 < lGeV 2 contribution to the leading hadronic contribution to the muon g — 2, and give 
these results in Table IIVI We will leave the discusion of these results to the next section. 

More interesting are the fits using the SXPT vector resonance expressions for the vacuum 
polarization. We have performed three different fits: 

Fit A: Only the tree-level vector contribution given in Eq. (1301 ). 




(37) 



m=0 
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Fit B: Eqs. (IT7|) and (!30|) which include the tree- level p and the one-loop pion and kaon 
terms, but not the one-loop terms in Eq. fl35|) . 



Fit C: The full one- loop expression, but absorbing the one- loop p analytic term into fy. 

Note that the last two fits give identical results for since the pion and kaon logarithms and 
the NLO analytic term just serve to re-scale the (tree level) value of fy as the momentum 
dependence is identical to the tree-level expression. 9 In other words, the overall factor 
multiplying the (p 2 + m 2 /) -1 term stays the same, only the value of fy changes. Were we 
to be more precise, we would calculate the complete one-loop corrections to fy up to and 
including 0(q A ) terms and include these in the fitting functions. As our goal is to extract 
a value for the muon g — 2, such a calculation, which for consistency also entails two loop 
diagrams in the pion sector, is beyond the scope of this paper. 

The results from these fits with the maximum value of q 2 = 1 GeV 2 are shown in Table HTT1 
The values for the decay constant change slightly from Fit A to Fit B when adding in the 
pions, and there is some tendency for the decay constant to decrease with quark mass. In 
other words, the one-loop corrections from the pion sector have only a small influence. Fit 
C shows the tree-level value of the decay constant changes dramatically due to the one-loop 
corrections in Eq. (|35|) . The fits are shown along with the data in Figure [TBI We take for 
the chiral scale A = 1 GeV in Fits B and C. 

The effect of decreasing the maximum value of q 2 in the fit to 0.5 GeV 2 is small. The 
central value of fy changes by less than one percent, about two percent, and less than one 
percent for mi = 0.0031, 0.0062, and 0.0124, respectively. Raising the maximum to 1.5 
GeV 2 induces larger changes, but the x 2 values of the fits are significantly worse. In all 



cases the vector masses are fixed to their nominal values from [32|. The polynomial fits are 
less stable, with very large errors for the smaller range, and significant undershooting of the 
data for the largest range. For these reasons we quote values of aj^ LO using the "best" fit 
range, < q 2 < 1.0 GeV 2 . 

The values for a^ LO resulting from the above fits are listed in Table IIVI and displayed 
in Figure [TH with only statistical errors shown. First, for the polynomial fits, we see a 
dramatic rise in a^ LO as we decrease the quark mass, and also as we increase the order of 



9 Note this would not be possible if the vector masses were not identical at this order, since then there 
would be a nontrivial q 2 dependence. 
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the polynomial; they are not stable in this sense. This is indicative of the calculation in 
general: the value of a^ LO is quite sensitive to the low momentum region, and hence the fit 
in this region, due to the nature of Eq. ([3]) and the smallness of the muon mass. The low q 2 
region is fit better as the order of the polynomial increases which increases the value of a^ LO 
but the errors on the fitted parameters increase such that the values for a^ LO also have large 
errors. Thus, it is preferable to use one of the physically motivated fitting functions. For 
Fits A and B, we note there is little difference in the final result for aj^ LO , as expected from 
the fit results themselves. As mentioned in the last section, there is no difference in a^ LO for 
Fits B and C since the one-loop corrections in Eq. 035p only rescale the tree-level value of 
fy. The statistical errors on a^ LO for the SXPT fits are much smaller than the polynomial 
ones, so we use the former fits from now on to quote our best values. 

From Figure [TBI we see that the fits tend to undershoot the lattice calculation of —U(q 2 ) 
for the lowest momenta for the smallest two quark masses, though within roughly a standard 
deviation. As for the polynomial fits, even small changes in the fits in this region lead to 
large changes in a^ L °. This undershooting behavior could represent real physics, or simply 
statistical and systematic errors. Certainly, the values of — II(g 2 ) at the smallest values of q 2 
are the most difficult to calculate, so the latter is more likely to be the case. The good fits 
obtained using chiral perturbation theory and the precisely measured meson masses from 
{32) ]. over a wide range of momentum, also back up this explanation. Still, the possibility 
that we have not accounted for an effect due to small quark mass and momentum remains, 
and must be further investigated in future calculations. In such calculations it is important 
to reduce the statistical error on these points as much as possible. One way to do this is to 
use a momentum source for each of the lowest momenta. This should have smaller errors 
than the point (-spilt) source used here, but is more complicated to implement and requires 
a separate propagator calculation for each momentum. Once the statistical errors on the 
very low q 2 region are reduced, one can begin to investigate systematics to tell whether the 
excess is an actual physical effect (which could increase the value of a^ LO significantly). 

In a similar vein, we should also check the numerics of our calculation. At the heart of the 
calculation is the quark propagator computation which is performed using the conjugate gra- 
dient algorithm to invert the lattice Dirac operator (matrix). The accuracy of the inversion 
is controlled by stopping the iterative solution found by the conjugate gradient algorithm 
when the norm of the so-called residual vector reaches some small size. The stopping crite- 
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ria used here are given in Table d and are comparable in size to those used in other current 
simulations. They are also much smaller than the criteria used to generate the ensembles 



32J. One test was run on a single 40 3 , mi = 0.0031, lattice to check whether the stopping 



criteria was adequate. We calculated the vacuum polarization using three stopping criteria, 
10~ 5 , 10~ 6 , and 10~ 7 , and compared the results. The difference in last two calculations was 
always less than 0.1 % except for two cases for the lowest momenta where the difference was 
still less than 1%, while the difference between the first two was always less than 1%, except 
for a handful of cases including one for the lowest momenta where the difference was neary 
10%. Thus we chose to use a stopping criteria of 10 -6 in this case (for all others we used 
10 -8 ), and we expect in all cases that the error in our results for the vacuum polarization due 
to a non-converged propagator solution is much smaller than the quoted statistical error. 

The other possible problem with the numerics is the numerical precision of the calcula- 
tion. In all cases the global sums in the quark propagator calculations were done in double 
precision and the rest of the calculation done in single (32-bit) precision, including the fourier 
transform of the vacuum polarization. This was also the procedure used in 32j , although in 
that calculation they did a comprehensive check for systematic errors arising from the use 
of single precision arithmetic (see Table VIII of 32J). Any differences due to precision were 
much smaller than quoted statistical errors, or even statistical errors that might be achieved 
by a "reasonable lengthening" of their run. 

Thus for now we attribute the excess of the vacuum polarization over the fit results at 
the lowest one or two momenta for the smallest two quark masses to statistics and possible 
systematics (but not numerics), and rely on the fits to chiral perturbation theory to provide 
an accurate determination of the anomalous magnetic moment. We also note that the chiral 
perturbation theory fits for all quark masses are very stable if we leave out the lowest two 
momenta, as would be expected from the relative size of the statistical errors on these points. 

We need to extrapolate the values of a^ LO at fixed quark mass to the physical point 
(which, in lattice units, is a light mass of about ami ~ 0.001 32|, |33|]). There are several 



is 



ways to do this. The simplest way is to extrapolate in the light quark mass. This i 
shown in Figure [151 with two different "phenomenological" fits. One is a linear fit in the 
lightest two masses, and the other is a quadratic fit to all three. The value from the linear 
fit is 

a HLO,,a<i = (713 ± i5) x iQ-io (38) 
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while the quadratic gives 

a HLo,^<i = ( 742 ± 21) x 10 -u ( 39 ) 

where the errors are statistical. Of course, there is no guarantee that (other) non-linear 
behavior does not set in as the two pion threshold is crossed, which we discuss below. 

The more correct method to reach the physical value is to fit all the lattice data simultane- 
ously, and extrapolate the vector decay constant and vector and pion masses to the physical 
light quark mass, with a single mass- and momentum- independent constant. The difficulty 
lies in the fact that such a fit requires computing the effects of the two pion threshold on 
the vector masses and decay constants, a clearly formidable task as the p meson is unstable 
for physical quark masses. In our study, on the 40 3 lattice with ami = 0.0031, the vector 
meson is just below the two pion threshold; however it shows clear non-linearity in the quark 
mass (see Figure [T6|) . Lacking knowledge of two pion effects in chiral perturbation theory 
to the specified order at which we are working, we could instead simply develop smooth 



parameterizations of the masses (for example, see [46j]) and decay constants and use them 
in our chiral perturbation theory formulae. But such a complicated procedure is not likely 
to yield more accurate results from this study than the simple extrapolation we have used. 

Using the continuum three-loop perturbation theory result matched to the ami = 0.0031 
lattice calculation at 2 GeV as described above, we find 

a HLO^>l = g x 10 -10 ( 40) 

for the high energy region, q 2 > 1 GeV 2 . There is some small error in this result made 
by running down to 1 GeV 2 instead of some higher value more compatible with the lattice 
calculation as seen in Figure [101 However, since the contribution to a^ LO is already so small 
starting at 1 GeV 2 and the difference between the perturbative curve and the lattice data 
is small even at 1 GeV 2 , we did not bother to fit the lattice data with a stepwise procedure 
up to such a high momentum. Note also that the perturbative curve is insensitive to the 
value of the light quark mass, so we did not attempt to extrapolate this contribution to the 
physical light quark mass. 

Adding the contributions from high and low momentum regions, the total value of the 
lowest order hadronic contribution to g— 2 is 

fl HLO = ( 721 ± 15) x 1Q -10 ( 41 ) 
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for the linear fit, and 



fl HLO = ( 74g ± 21 ^ x m -io^ ^ 



for the quadratic. Our result is larger than the currently accepted theoretical numbers [3|, |4j, 
692.4 (5.9) (2.4) xl(T 10 and 711.0 (5.0) (0.8) (2.8) xlO" 10 , which use only e + e~ data and 
include r decay data, respectively. Considering our statistical errors and the systematic 
uncertainties described below, the calculations should be viewed as compatible for now. 



VI. DISCUSSION 

This calculation of aj^ LO has focused on a detailed study of precisely fitting the lattice 
computed vacuum polarization to obtain the value of g — 2 for several fixed quark masses. 
This has up to now only been done at a single lattice spacing, and as such, there are still 
discretization effects that could arise that we have not taken into account. Simulations 
using the coarser MILC lattices are underway j^, and new lattices with a smaller lattice 
spacing than the 0.086 fm used here are being generated by the MILC collaboration, so the 
continuum limit of the vacuum polarization can be taken. One way to gauge the size of the 
non-zero lattice spacing errors in our calculation, without actually repeating the calculation 
at another lattice spacing, is to use the naive continuum momenta 2irn/L everywhere in the 
analysis instead of the sin function (2sin7m/L appears in Ward-Takahashi identity so is the 
correct form to use). The resulting difference in a^ LO is within the statistical variation of 
the previous result. Other evidence that non-zero lattice spacing errors are under control is 
the agreement between the lattice calculation of IT(g 2 ) and continuum perturbation theory, 
as already discussed in the previous section. 

We also have not explicitly taken into account possible finite volume effects. The quark 
masses in these configurations are small enough to have small but not negligible finite volume 



corrections 32j, |33(. Most likely, however, these corrections are small compared with the 
statistical errors that exist in the data. In terms of our model used to fit the data, the 
dominant contribution is from a tree-level photon-p-photon diagram, where there will be no 
corrections coming from finite volume (other than the finite volume correction to the vector 
mass itself which is expected to be quite small). Finite volume corrections first appear in 
the one-loop diagrams and thus will be suppressed relative to the dominant contribution. 
Another source of systematic error lies in the mismatch of the fixed value of the strange 
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quark mass on these lattices and the physical strange quark mass which is actually about 10 
percent smaller. Although we have assumed that the effect of this mismatch is negligible, this 
may not be the case. We can see from the lattice data that there is significant quark mass 
dependence for the vacuum polarization for low values of momenta. Note that correcting 
for this effect will tend to increase the value of a^ L °. 

Finally, as discussed already, the extrapolation to the physical light quark mass is not 
fully under control. Even though the extrapolations shown in Figure [15] appear sensible 
and are not large in the quark mass, the physical vector particles are unstable with non- 
zero widths, so these effects should be properly included in our analysis of resonance chiral 
perturbation theory. Such a calculation was beyond the scope of the present study as a 
consistent treatment of the pion and vector sectors entails higher order calculations. 

In this context it interesting to discuss the effects due to quenching since it is in this 
approximation that the vector particles are stable. In Figure [17] we show the quenched 



value of — II(g 2 ) computed on an ensemble of lattices |32| with the same lattice spacing 
and volume as the ami — 0.0062, 2+1 flavor ensemble. The quenched and 2+1 flavor 
results agree at large q 2 , but a noticeable deficit occurs for the quenched case as q 2 — > 0. 



The resulting value o: 
quenched calculations 



a^ LO is significantly lower (see Table HVl) . consistent with previous 
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151 ] . Presumably this effect is caused by the sea quarks. 
Even though the p mass is always below the two pion threshold in our calculation 10 , 
these two pion states contribute to the vacuum polarization, albeit with larger invariant 
mass, and on the lattice, a discrete spectrum. In our SXPT analysis two pion effects come 
only from the direct coupling to the photon; the two pion coupling to vector mesons is 0(p 4 ) 
and was not included. Because this coupling is higher order, we also expect that the effect 
of the non-zero width of the vector mesons on aj^ LO is suppressed. The difference in the 
two cases is reflected, in the fits, through the vector mass itself which is significantly larger 



in the quenched case 32] 11 and the decay constant which is also larger (see Table [ITT])- A 
more complete understanding of these effects is wanting, and must await both simulations at 
smaller quark masses where the vectors are unstable and higher order calculations in chiral 
perturbation theory, some of which have been done already in continuum chiral perturbation 



in 



Th e pi ons are required to have relative momentum, so the lowest pion energy is roughly yl (2irn/L) 2 + 



32j | only light quark masses 0.015 and 0.03 were computed in the quenched case, so we simply extrap- 



In 

olated linearly in am; to obtain the value of the vector mass at 0.0062. 
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theory 431 ] . We conclude this dicussion by emphasizing that there is a significant effect, 
clearly visible in the lattice calculations (see Figure H7|) . from unquenching which tends to 
increase the value of a^ L °. 

H 1 



VII. CONCLUSION 

In this study we presented a value of the lowest order hadronic contribution to the muon 
anomalous magnetic moment computed in 2+1 flavor QCD. The statistical errors on our 
results at fixed light quark mass are small, possibly underestimated, and comparable with 
the errors on the accepted value of these contributions computed using the experimentally 
measured total cross sections for e + e~ annihilation to hadrons and tau lepton decay to 

n 

hadrons[3|. The central value found by extrapolation to the physical light quark mass is 
slightly larger than this value, with statistical errors about 2-3 times larger. Within the 
systematic uncertainties in our calculation which were discussed in detail in the previous 
sections, the values should be seen as compatible. 

First, there is some uncertainty in our fitting procedure. Polynomial fits require a high 
degree to fit the low momentum region of the vacuum polarization so that the statistical 
errors on the parameters become large. The results are also somewhat unstable to changes 
in the fit (momentum) range. 

A more precise and stable fitting procedure was found by appealing to chiral perturbation 
theory, including lattice spacing effects. To our knowledge, this is the first time the hadronic 
vacuum polarization computed in chiral perturbation theory has been given in the literature. 
While the (staggered) chiral perturbation theory formula depends only on one (unknown) 
low energy constant and thus offers a useful fit function, we found that it does not represent 
the lattice calculation; to fit the data well requires the inclusion of the vector particles as 
well, through resonance chiral perturbation theory. This was already indicated by quenched 
calculations. The additional free parameter is the vector decay constant which shows modest 
quark mass dependence at tree level and large effects due to one-loop pion and kaon tadpole 
graphs. The coupling of the vector meson to two pions which contributes to both mass 
renormalization and the non-zero width of the p meson gives rise to a one-loop graph which 
is however 0(p A ), and hence was not included in this work since a consistent treatment in 
the pion sector would require a two-loop computation. We also emphasize that a direct 
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comparison to a "matched" quenched calculation showed a significant increase in the 2+1 
flavor vacuum polarization at low momenta, presumably due to sea quark effects. 

The vacuum polarization computed on the lattice was seen to be compatible with con- 
tinuum weak-coupling perturbation theory over a wide range of momenta, between roughly 
3 and 8 GeV 2 . This indicates that lattice spacing errors may be under control in this study 

We have used the fits of the vacuum polarization to chiral perturbation theory to quote 
our final results. Remaining fitting uncertainties can be reduced in future calculations by 
reducing the statistical errors on the smallest momentum values of IT(g 2 ) and by including the 
0(p 4 ) terms in the chiral perturbation theory analysis. The size of the statistical error can 
be further reduced by increasing the ensemble size, computing with momentum sources, and 
computing at smaller valence quark mass. These require modest increases of computational 
resources compared to those used here. Reducing the sea quark mass is obviously very costly, 
but not unforeseeable. 

There are additional, as yet unquantified, systematic errors in our calculation, notably 
effects due to omission of disconnected quark diagrams giving rise to flavor singlet contribu- 
tions, non-zero lattice spacing and unphysical quark masses. The latter were discussed at 
some length in the last section. The main point is that the two pion threshold may cause 
non-linearities in the light quark mass extrapolation which can be addressed by simulating 
below the threshold or, at least, by including higher order terms in the chiral expansion. It is 
also crucial to take the continuum limit and compute the vacuum polarization with different 
lattice fermions in light of the current understanding of the non-local nature of staggered 
fermions at non-zero lattice spacing. While these effects will be addressed in future calcu- 
lations, the current study of the muon g — 2 on the lattice in 2+1 flavor QCD appears quite 
encouraging. 
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TABLE I: Parameters for the MILC lattices and propagator calculations in this study, "e.g. stop 
res." is the stopping criterion, used for the residual vector in the conjugate gradient subroutine 
of the MILC code, for the valence quark propagator. The stopping criteria used to generate the 
lattices are given in Table I of 



32]. In the "confs." column total number of configurations are given 
for both propagator source time slices, t = and 48. (3 = W/g 2 is the bare gauge coupling for the 



tadpole-improved Symanzik gauge action used to generate the 
and vector meson masses used in this study were taken from 



attice ensembles 



3J] and 



3Jj. The vr, K, 



(ami = 0.0031). 



a (fm) size f3 ami ar "s ar ™val e.g. stop res. confs. 



am, 



amj< 



amy 



0.086(2 
0.086(2 
0.086(2 
0.086(2 
0.086(2 
0.086(2 
0.086(2 
0.086(2 



28 3 x 96 8.40 - - 0.031 
28 3 x 96 8.40 - - 0.0062 
28 3 x 96 7.11 0.0124 0.03 1 0.031 
28 3 x 96 7.11 0.0124 0.03 1 0.0124 
28 3 x 96 7.09 0.0062 0.03 1 0.031 
28 3 x 96 7.09 0.0062 0.03 1 0.0062 
40 3 x 96 7.08 0.0031 0.03 1 0.031 
40 3 x 96 7.08 0.003 1 0.03 1 0.0031 



10~ 8 
10~ 8 

io- 8 

10~ 8 

io- 8 
io- 8 
io- 8 

10" 6 



84/54 
111/54 
129/129 

129/129 0.20638 (18) 0.27209 (18) 0.4173 (13) 
42 

249 0.14794 (19) 0.25319 (19) 0.3895 (28) 
51/51 

161/153 0.10521 (8) 0.24175 (15) 0.357 (5) 
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TABLE II: Fit parameters for the cubic and quartic polynomials for the 2+1 flavor value of — Il(q 2 ). 
The fit range was taken to be < q 2 < 1 GeV 2 in each case. The jackknife estimates of the errors 
are statistical only, and the x 2 values are from uncorrelated fits. 



order 


ami b 


bi 


b 2 


b 3 


bi X 2 /dof 


3 


0.0124 0.0972(14) 


-0.051(7) 


0.038(13) 


-0.014(7) 


3.5/16 


4 


0.0124 0.0978(20) 


-0.058(14) 


0.062(39) 


-0.049(45) 


0.017(19) 3.4/15 


3 


0.0062 0.1025 (14) 


-0.0615 (63) 


0.045 (10) 


-0.015(5) 


19.3/7 


4 


0.0062 0.1058 (27) 


-0.092 (19) 


0.137(48) 


-0.124(51) 


0.045(19) 16.3/16 


3 


0.0031 0.1075 (9) 


-0.0754 (35) 0.0636 (54) 


-0.0241 (27) 


25.6/35 


4 


0.0031 0.1102(13) 


-0.1044(89) 


0.157(23) 


-0.142(26) 


0.051(10) 17.4/34 



TABLE III: Fit parameters for SXPT formulae for the 2+1 flavor value of — U(q 2 ). The first 
row is for the quenched case discussed in the text. The fit range was taken to be < q 2 < 1 
GeV 2 . The jackknife estimates of the errors are statistical only, and the value of x 2 /dof is from 
an uncorrelated fit. The meson masses were fixed to the values given in Table HI The staggered 
meson mass splittings used in Fit B and Fit C (not shown) are found in 321 ] . 



am x f™ A (MeV) A FitA x 2 /dof /£ itB (MeV) A FitB x 2 /dof ff tC (MeV) A FitC x 2 /dof 
0.0062 209.9 (2.0) 0.0410 (6) 23/18 

0.0124 192.8 (1.8) 0.0445 (6) 4.4/18 188.6 (2.0) 0.0421 (7) 4.0/18 117.2 (1.2) 0.0420 (6) 4.0/18 
0.0062 186.8 (1.7) 0.0453 (5) 18/19 181.6 (1.8) 0.0422 (5) 18 /19 115.7 (1.0) 0.0422 (5) 18/19 
0.0031 175.4 (1.1) 0.0474 (3) 28/37 169.9 (1.1) 0.0436 (3) 25/37 111.5 (7) 0.0436 (3) 25/37 
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TABLE IV: Results for a^ LO x 10 for the various fits described in the text. Errors are jackknife 
estimates and statistical only. The quenched results correspond to light valence quark mass 0.0062 
and strange valence quark mass 0.031. 



Fit quenched 


ami = 0.0124 


ami = 0.0062 


ami = 0.0031 


Poly 3 381 (63) 


370(49) 


445(43) 


542(24) 


Poly 4 588 (142) 


410(91) 


639(123) 


729(59) 


A 366.6 (7.0) 


412.3 (7.8) 


516.0 (9.5) 


646.9 (8.1) 


B 


403.9 (7.8) 


502.1 (9.5) 


628.0 (8.1) 


C 


403.9 (7.8) 


502.1 (9.5) 


628.0 (8.1) 
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FIG. 1: The lowest order diagram giving rise to QCD corrections to the muon g— 2, at 0(a 2 ). The 
blob represents all possible hadronic states. 




FIG. 2: The two Wick contractions for the vacuum polarization. The dots are insertions of the 
vector current J^x). Only the first has been calculated in this study. 
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FIG. 3: The two one-loop diagrams contributing to the photon vacuum polarization in SXPT. The 
solid line can be either a charged pion or a charged kaon. 




FIG. 4: Quark-level vertex coming from the one-photon term in the Lagrangian. The solid lines 
are quarks. 




FIG. 5: Quark-level vertex coming from the two-photon term in the Lagrangian. As in Figured! 
the solid lines are quarks. 
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FIG. 6: Comparison of Eq. (|17p to the lattice vacuum polarization, where we have matched the two 
at the smallest value of q 2 to fix the overall constant. The open circles are the datapoints (errors 
not shown) for the lightest quark mass used, ami = 0.0031. The solid line is the SXPT curve. 
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FIG. 7: Tree-level diagram for a vector meson which contributes to the vacuum polarization. The 
double- line can be either the p° or the cog. 



FIG. 8: The one-loop diagram involving the vector mesons. The solid line is either a pion or kaon, 
and the double solid line is the vector meson. 
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q 2 (GeV 2 ) 

FIG. 9: Minus the vacuum polarization for 2+1 flavors of quarks for each light quark mass studied 
in this work, 0.0031 (diamonds), 0.0062 (squares), and 0.0124 (circles). The insert shows a blow 
up of the important low q 2 regime. The strange quark mass is fixed to 0.031 in each case. 
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FIG. 10: Comparison of the 2+1 flavor vacuum polarization computed using the lattice with 3-loop 



continuum perturbation theory in the MS scheme 



441 ] . The solid line is forced to match the lattice 



calculation at 2 GeV through a simple additive shift, mi = 0.0031 and m s = 0.031. The quark 



masses have been converted to the MS scheme using the matching factor in 



3J 
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FIG. 11: Same as Figure [10] but for mi = 0.0124. 
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FIG. 12: Cubic (dashed) and quartic (solid) fits to -II(g 2 ) for ami = 0.0031 (diamonds), 0.0062 
(squares), and 0.0124 (circles). The strange quark mass is fixed to 0.031 in each case. 
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FIG. 13: SXPT fits to U(f) for the three light masses, ami = 0.0031 (diamonds), 0.0062 (squares), 
and 0.0124 (circles). The strange quark mass is fixed to 0.031 in each case. The solid lines 
correspond to Fit B, the dashed to Fit A, as described in the text. 
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FIG. 14: The hadronic contribution to the muon anomalous magnetic moment from polynomial 
and SXPT fits to the vacuum polarization. The statistical uncertainties for the cubic (squares) and 
quartic (circle) fits are much larger than the SXPT one. 
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FIG. 15: Linear and quadratic extrapolations (bursts) of the hadronic contribution to the muon 
anomalous magnetic moment to the physical point (ami ~ 0.001). The left triangle is the currently 



accepted value, computed from the experimental cross-section for e + e — > hadrons 



IB. 
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FIG. 16: The vector meson mass (GeV). Non-linearity is visible for the lightest quark mass studied 
here, and perhaps even the second lightest mass. The horizontal line depicts the two pion energy 
for the largest lattice and smallest quark mass in our study. The mass values are from the MILC 
collaboration, 



3j and 



4$ (ami = 0.0031). 



13 



0.12 



~i 1 1 1 1 1 1 1 1 1 1 1 1 r 

0.12 



t 1 1 1 1 1 r 



0.08 



0.06 



0.04 



~\ — i — i — i — I — i — i — i — r 




■oo-oq., 



■<5>.000 



J I I I I I I I L 



J I I I I I I I L 



q 2 (GeV 2 ) 



FIG. 17: The vacuum polarization computed on the quenched ensemble (circles) compared to the 
2+1 flavor result. The light quark mass in each case is 0.0062, and the quenched gauge coupling 
was tuned to match lattice spacings. The solid lines are fits to Eq. (|30p . 
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